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Abstract. A large population of fragile, wide (> 10 3 AU) binary systems exists in 
the Galactic field and halo. These wide binary stars cannot be primordial because of the 
high stellar density in star forming regions, while formation by capture in the Galactic 
field is highly improbable. We propose t hat these binary systems w ere formed during 
the dissolution phase of star clusters (see iKouwenhoven et al.ll20lol for details). Stars 
escaping from a dissolving star cluster can have very similar velocities, which can lead 
to the formation of a wide binary systems. We carry out A^-body simulations to test 
this hypothesis. The results indicate that this mechanism explains the origin of wide 
binary systems in the Galaxy. The resulting wide binary fraction and semi-major axis 
distribution depend on the initial conditions of the dissolving star cluster, while the 
distributions in eccentricity and mass ratio are universal. Finally, since most stars are 
formed in (relatively tight) primordial binaries, we predict that a large fraction of the 
wide "binary stars" are in fact higher-order multiple systems. 



1. Wide binary systems in the Galactic field and halo 

The large m ajority of stars are thought to form as part of a binary or multiple stellar 



system (e.g ., Duquennoy & Ma yor! 1 19911 : iFischer & Marcylll992l : iGoodwin & Kroupa 
20051 : IKouwenhoven et al.ll2005L l2007b~ The general consensus is that most sta r form 



in embedded star clusters and loosely-bound associations (e.g., lLada & L ada 2003 



Bastian 201 1), which initially exhibit a significant amount of substructure (e.g jAllison et al. 



2009). Following proto-star formation, the properties of the binary population ev olve 
over time, primarily due to the effects pre-ma in sequence evolution (IKroupall 19951) an d 
dynamical interactions with other stars (e.g., iHeggie & Hutll2003l: iMarks et all 1201 lh. 



Most star clusters dissolve within 10-50 Myr after their formation (see lde Grijs & Parmentier 
2007, and references therein). The field star population is theref ore thought to b e the 
result of a mixture of stars originating from different star clusters dGoodwinll2010l) . 
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Over the last decades a significant number of wide (> 10 3 AU) binary systems have 
been d iscovered (see Fig.Q]). In the log-normal period distribution of lDuquennoy & Mayor! 
dl99ll) . for example, ~ 15% of the binary systems have a semi-major axis larger than 
10 3 AU. Individual wide binary systems are often identified in proper motion studies, 
and occasionally combined with parallaxes, radial velocity measurements and back- 
groun d star statistics (e.g., Mak arov et al.ll2008l : lOuinn & Smithll2009l : IShaya & Oiling! 
201 1| . and numerous others). The o verall properties of the wide b inary population can 



also be obtained statistically (e.g., lLonghitano & BinggelJ |2010|) . Wide binary sys- 
tems are extremely fragile, and those wider than 0.1 - 0.2 pc are easily destroyed in 
the Galactic field (see Fig. [j}. This upper limit can be explain ed by interactions with 
other stars, molecular c louds, and the Galactic tidal field (e.g., Rett erer & Kindl 19821 : 



Jiang & Trema ine 2010). The properties of wide systems in the Galactic field are also 



used to constrain the properties of hypoth esized dark components (e.g.. lOuinn et al. 



used to constrain tne properti es ot nyp otnesized 
2009; Al len et al.ll2007l : iHernandez & Led l2008\ 



Wide binary system^] cannot have formed as primordial binaries in star clus- 
ters, simply because their orbital separation is comparable to the size of a typical 
embedded cluster. Moreove r, the typical size of a star forming core is ~ 10 4 AU 



(IWard-Thompson et all 120071) . which sets an absolute maximum to the size of a pri- 
mordial wide binary system. Even if they were somehow able to form, they wou ld be 
destroyed immediately due to stellar encounters dKroupall200lUParker et al. [2009). 



The fact that it is not possible to form primordial binary systems with semi-major 
axes in the range 10 3 AU-0.1 pc implies that wide binaries are formed at a later stage, 
as a result of dynamical interactions between stars. Energy conservation implies that 
two stars on an initially unbound orbit will remain unbound. Capture is therefore only 
possible when energy is removed from the system, for example by a third star that is 
present during the encounter. The formation rate Nb of binary systems via three-body 
encounters is given by 

G 5 M 5 n 3 

N B = 0.75 5— , (1) 



(IGoodman & Hutl ll993). where G is the gravitational constant, M is the typical stellar 
mass, n the stellar number density, and cr the velocity dispersion. The value of Nb is 
negligible for stars in the Galactic field and halo. On the other hand, capture is possible 
in the dense cores of star clusters, but this will never result in the formation of long-lived 
wide bin ary systems, due to the cro wded stellar environment. 

In [Kouwenhoven et all (1201 Oh we proposed that wide binary systems form during 
the dissolution phase of star clusters . This mechanism can result in a significant pop - 



ulation of binary systems (see also Moecke f & Batdl201dl : iMoeckel & Clarkel 12011 ). 



In this scenario, an unbound pair of escaping stars can form a binary system when 
their relative velocity is smalfl. Our Af-body simulations (see below) result in a pop- 
ulation of wide binary systems with semi-major axes comparable to their initial sepa- 
ration, a thermal eccentricity distribution, and a mass ratio distribution resulting from 
(gravitationally-focused) random pairing of components from the initial mass function. 



'Following Kouwenho ven et all d2010h . we define systems with an orbital separation in the range 10 3 AU 
< a < 0.1 pc as wide binary systems (or wide multiple systems). 

2 A large populati on of comets may also be captured b y stars via a similar mechanism, during cluster 
dissolution (e.g., Eg gers et alll997l : lLevison et alj20ioh 
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Figure 1. The observed semi- major axis distribution for wide binary systems , 
compiled from the catalogues of [ Duquennoy & Mayor ( 1991) :IClose et al.| (| 1990): 
iLepine & Bongiornol j2007l) : IChaname & Gouldl d2004l) and lPoveda et alj d2007l) . 



2. Method, initial conditions and Af-body simulations 

To test our hypothesis that wide binary systems are formed during the dissolution 
phase of star clusters, we ca rry out A-body simulations using the STARLAB pa ckage 



( Portegies Zwart et al.ll200lh . For each star cluster were draw N stars from a liCroupa 
(2001) mass distribution in the range 0.1 -50M Q . We study the properties of the result- 
ing binary population as a function of the number of member stars N (10 <N< 1000), 
the size R of the star clusters (0.1 pc < R < 1 pc), and the primordial binary fraction 
B (0% < B < 100%). We additionally vary the virial ratio Q = E K /E p of the clus- 
ter, where and E p are the kinetic and potential energy of the cluster, respectively. 
We study the cases 2 = 1/2 (cluster in virial equilibrium) and 2 = 3/2 (expanding 
star cluster). We carry out simula t ions f or two different stellar density distributions 
of the stars: (i) spherical iPlummerl (Il91 ll) models, a nd (ii) models with substructu red 
initial conditions with fractal parameter a = 1 .5 (see iGoodwin & Whitworthll2004l for 
details). 

Simulations are carried out until the modeled cluster has completely dissolved. At 
the end of each simulation we determine the properties of the resulting binary popu- 
lation. We identify a pair of stars as a binary system when (i) their binding energy is 
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Q = 1/2 
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Figure 2. The semi-major axis distribution (left), the correlation between mass 
ratio q and semi-major axis a (middle) and between primary mass M\ and mass ratio 
q (right). The properties of the orbits of binary systems and higher-order multiple 
systems are indicated with the dots and triangles, respectively. For each multiple 
system with n stellar components, we have included all n — 1 orbits. Results are 
shown for 50 Plummer models with N = 1000 and R = 0.1 pc, and virial ratios of 
2=1/2 (top) and 2 = 3/2 (bottom). The vertical dashed lines indicate a = 10 3 AU 
and a = 0.1 pc, respectively. The dashed curve in the right-hand panel indicates the 
minimum mass ratio <7 m i n (Mi) = M m ; n /Mi. 



negative, and (ii) both stars are each others mutual nearest neighbor. In our analysis 
we only consider binary systems with a semi-major axis a < 0.1 pc, as wider pairs are 
rather easily destroyed due to stellar encount ers. We also id e ntify hierarchical multi- 
ple stellar systems (> 3 stars) and impose the IValtonen et all d2008h stability criterion 
^outMn > Q st for multiple stellar systems, where a m and a out are the inner and outer 
orbits of a (sub)system, respectively, and Q st « 3 - 10 is a stability parameter which 
depends on the orbital configuration. For wide higher-order system with (outer) orbital 
periods of order ~ 1 Myr, this corresponds to a stability timescale of several billion 
years. 



3. The wide binary fraction, orbital elements and higher-order multiplicity 

A summary of the results of our Af-b ody simulations is listed be low. For an extensive 
description of the results we refer to iKouwenhoven et al. 1 fcoich . To illustrate the re- 
sults, we show the orbital properties of a selected sample of the simulations in Figs. |2] 
and [3] 



(a) The wide binary fraction. After cluster dissolution, the wide binary fraction (i.e., 
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Figure 3. Same as Fig. [2] but now for substructured models (fractal parameter 
a = 1.5) with N = 1000 and R = 0.1 pc , and virial ratios of Q = 1/2 (top) and 
2 = 3/2 (bottom); in each case fifty realisations have been simulated. 



the fraction of wide binary systems as compared to the total number of systems) ranges 
between 1% and 30% for an individual star cluster. The exact value depends the prop- 
erties of the star cluster at the time of dissolution. The structure of the star cluster at the 
moment of dissolution, as well as the number of stars, affects the final number of wide 
binary stars. Substructured star clusters (e.g., Fig. [3]) generate significantly more wide 
binary systems than spherical star clusters (e.g., Fig. O. The wide binary fraction in- 
creases with decreasing cluster membership and with increasing initial virial ratio. The 
wide binary population in the Galactic field (~ 15%) results from a mixture of wide 
binary systems formed from different types of star clusters. Its properties can therefore, 
in principle, be used to constrain the properties of young star clusters. 

(b) The semi-major axis distribution. The resulting semi-major axis distribution for 
wide binary systems is mainly in the range (0.1 - l)R, where R is size of the star cluster 
at the moment of dissolution. The semi-major axis distribution of the newly formed 
binaries typically shows two peaks: a dynamical peak at small values of a, resulting 
from three-body interactions, and an dissolution peak of wide binary systems formed 
during the dissolution phase of the star cluster (this is clearly shown in the left-hand 
panels in Fig. |3]>. The ratio of the number of binary stars in the dynamical peak and the 
dissolution peak depends on the initial conditions of the star cluster (see above). 

(c) The eccentricity distribution. The capture process which results in wide binary 
formation is chaotic. The eccentricity distribution for wide bin ary systems is therefore 
expected to be thermal: f(e) = 2e for < e < 1 (He ggie|[i975l) . which is confirmed by 
the simulations. 
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Figure 4. Most stars form as a member of a primordial binary system. Many wide 
"binary" systems formed during the star cluster dissolution process are therefore 
expected to be higher-order multiple systems. The relative multiplicities of wide 
systems can in principle be used to constrain the binary fraction Bq at the time of star 
cluster dissolution. 



(d) The mass ratio distribution. The mass ratio distribut ion for wide binary system s re- 
sults from gravitationally-focused random pairing (e.g.. lKouwenhoven et alJl2009L and 
references therein) of the individual components. This implies that wide binaries with 
a high-mass primary star have a small mass ratio, while those with a low-mass binary 
have a high mass ratio (see Figs. |2] and [3]). In addition, the wide binary fraction slowly 
increases with increasing primary star mass. 



(e) Orientation of the orbits. The orientation of the stellar spins of the two stars in 
a wide binary system are randomly aligned. In the case of a wide multiple system, the 
orbital orientations of the inner orbits are also randomly aligned with respect to each 
other, and with respect to the orbit of the wide orbit. 



if) Implications for higher-order multiplicity. A significant fraction of star form as 
part of a primordial binary system. Both components of a wide "binary" system are 
therefore expected to be binary themselves. Recent observation s suggest indeed that 
wide "binary" system s are frequently triple or quadruple sys tems dMakarov et al.ll2008l : 
Mamajek et all l201fi iFahertv et aLll2010l : lLaw et all |2010|) . The multiplicity ratios 
among wide systems can therefore be used to constrain the primordial binary fraction, 
or more specifically, the binary fraction at the moment of star cluster dissolution (see 
Fig. 11. 



There is an ongoing debate about the wide binary formation mechanism itself. 
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(2010) shows that a pair of (previously unbound ) stars can form a wide b inary system 
during the dissolution phase, while the study of iMoeckel & Clarke] (1201 lh shows that 
a small (but transient) population of wide binary systems is always present in a star 
cluster, and that this wide binary population is frozen in when a star cluster dissolves. 
It may well be possible that both mechanisms contribute to the formation of the wide 
binary population in the Galactic field and halo. 

4. Summary 

Approximately 15% of the known binary systems in the Galaxy have an orbital sepa- 
ration larger than 10 AU. These systems cannot be primordial, simply because their 
orbital separations are comparable to the size of young embedded clusters. Moreover, 
if they were able to form in such environments, they would immediately be destroyed 
by dynamical interactions with other stars. Dynamical capture in the Galactic field or 
halo is highly improbable due to the low stellar density and high velocity dispersion, 
and cannot explain the observed wide binary population either. 

We pro pose that wide binary sys tems form during the dissolution phase of star 
clusters (see lKouwenhoven et alj|2010l for details). In this scenario, an escaping pair of 
stars with a small relative velocity can form a wide binary system. TY-body simulations 
confirm this hypothesis, and allow us to predict the prevalence and orbital properties of 
the wide binary population (§ [3]>, and the fraction of triple and quadruple stars among 
wide systems (see Fig.@]). These predictions can be tested observationally, in particular 
those for the mass ratio distribution and the higher-order multiplicity. 
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